function [intZ,intZsquare]=intOU(t,T,Npath)



w = randn(5000,Npath);
getOU = oupath(w,1,1,1,1);

intZ = (sum(getOU,1)-getOU(1,:)/2-getOU(5001,:)/2)*(T-t)/5000;
getOU = getOU.^2;
intZsquare = (sum(getOU,1)-getOU(1,:)/2-getOU(5001,:)/2)*(T-t)/5000;